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ABSTRACT 

In this paper, we examined XMM Newton EPIC spectra of 14 ultra-luminous X-ray sources (ULXs) 
in addition to the XMM RGS spectra of two sources (Holmberg II X-1 and Holmberg IX X-1). We 
determined oxygen and iron abundances of the host galaxy's interstellar medium (ISM) using K-shell 
(O) and L-shell (Fe) X-ray photo-ionization edges towards these ULXs. We found that the oxygen 
abundances closely matched recent solar abundances for all of our sources, implying that ULXs live 
in similar local environments despite the wide range of galaxy host properties. Further, the ISM 
elemental abundances of the host galaxies, as indicated from the 0/H values obtained from the X-ray 
spectral fits, are in good agreement with the 0/H values obtained by Garnett (2002) from studies 
of H II regions in the same external galaxies. We find roughly solar 0/H values independent of the 
host galaxy luminosity, star formation rate, or the ULX location in the galaxy. Also, we compare the 
X-ray hydrogen column densities (nj^) for 8 ULX sources with column densities obtained from radio 
H I observations. The X-ray model t\h values are in good agreement with the H I values, implying 
that the hydrogen absorption towards the ULXs is not local to the source (with the exception of the 
source M81 XMMl). In order to obtain the column density and abundance values, we fit the X-ray 
spectra of the ULXs with a combined power law and one of several accretion disk models. We tested 
the abundances obtained from the XSPEC models bbody, diskbb, grad, and diskpn along with a 
power law, finding that the abundances were independent of the thermal model used. We comment 
on the physical implications of these differe nt model fits. We also note t hat v ery deep observations 
allow a breaking of the degeneracy noted bvl Stpbbart. R.oberts. fc Wilmsl |)200(tI) favoring a high mass 
solution for the absorbed grad + power law model. 
Subject headings: ISM: general — ISM: abundances — X-rays:ISM 



1. INTRODUCTION 

Long exposure XMM-Newton observations of nearby 
galaxies offer new opportunities to study various proper- 
ties of ultraluminous X-ray source (ULX) spectra. ULXs 
are bright, non-nuclear X-ray sources with X-ray lu- 
minosities > 3 x lO'^^ergs"^. In a previous paper 
(jWinter. Mushotzkv. fc Revnoldsl l2005|) (Paper 1) , we 
analyzed ULX spectra from 32 nearby galaxies. Based on 
spectral form, luminosity, and location within the optical 
host galaxy, we classified a population of high/soft state 
and low/hard state ULXs. For the ULX sources with the 
greatest number of counts, the high signal-to-noise XMM 
observations can be analyzed with more than the simple 
schematic models of our first study (an absorbed power 
law and blackbody model for the high-state and an ab- 
sorbed power law model for the low state). Particularly, 
the spectra can be used to investigate the properties of 
absorption along the line of sight to ULX sources. 

Many similar studies have been done within our 
own Milky Way, where X-ray absorption models have 
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been used to determine column densities and abun- 
dances of the interstellar medium (ISM). The proce- 
dure involves spectral fits of absorption features in 
bright, background X-ray sources. Successful determi- 
nati ons have been made using backgro und galaxy clus- 
ters |I3aumgartjier_&_^£ushotzkv"2006') and X-ray bina- 
ries |luctt. Sc hulz. fc Ch akrabartv 2004). These stud- 
ies made use of a bright, X-ray source as a background 
through which they can observe the 542 eV absorption 
edge produced by photo-ionization of the inner K-shell 
electrons of oxyg en. Analogous studies have been used 
in the radio fsee iDickev fc Lock man (1993)) optical 
regime, using quasars, supernovae, or stars as a back- 
ground for hydrogen absorption and 21-cm emission, as 
a means to measure hydrogen column densities and metal 
abundances. 

In this study, we extend the X-ray absorption studies 
to external galaxies using ultra luminous X-ray sources. 
Due to their extreme brightness in the X-ray regime 
and their non-nuclear location in external galaxies, these 
sources are ideal for probing the ISM of their host galax- 
ies. Typical ULXs, from our lPaoer ll studv. have Galac- 
tic line-of-sight column densities of a few 10^" cm~^ 
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ijDickev fc LockmaDlll99d(l and measured X-ray colu mn 
densities greater than 10 cm" -2 (Fig. 8 of lPa,r)er ID for 
the combined ULX and host galaxy. Thus, if the local en- 
vironment of the ULX contributes little absorption, the 
X-ray column density is dominated by the host galaxy. 
One goal of this study is to determine whether this ab- 
sorption is that of the host galaxy or local ULX envi- 
ronment. Therefore, we compare the X-ray measured 
hydrogen column density with H I measurements from 
alternate methods. 

In addition to the brightness of ULXs and the rel- 
atively small Milky Way contribution to their X-ray 
hydrogen column densities, their well characterized X- 
ray spectra make ULXs ideal for measuring absorp- 
tion features of the ISM. Bright ULXs (e.g. NGC 
1313 X-1) typically have spectra that are well-fit by 
an absorbed multi-component blackbody and power law 
model. However, there is discussion over whether this 
standard model is the most physical model for the UL Xs 
(see, fo r example, Stobbart. Roberts. & Wilms (2006^ or 
IGoncalves &: Soria (2006)). Different models applied to 
the base ULX spectra can affect the absorption measure- 
ments, particularly in the softer part of the spectrum. 
Thus, in this paper we investigate the effect different soft 
component models have on the X-ray measured hydro- 
gen column density and elemental abundances (through 
the oxygen K-shell edge at 542 eV and the iron L-shell 
edge at 851 eV). 

We use high signal-to-noise XMM Newton observations 
of ULXs to measure hydrogen column densities and ele- 
mental abundances of oxygen and iron. Located in ex- 
ternal galaxies, the X-ray spectral resolution of avail- 
able ULX spectra is not as good as those of Galactic 
X-ray bin aries, which often have grating spec tra avail- 
able fe.g. lJuett. Schulz. fc Chakrabartvl l)2004() 'l. There- 
fore, in order to be able to distinguish the oxygen K- 
shell edge as well as the iron L-shell edge located at 
851 eV, we needed observations with a large number of 
counts (« 5000 counts). The X-ray observatory XMM- 
Newton, having a larger collecting area than Chandra, 
provides the counts necessary in order to conduct this 
study. Further, with recent 100 ks XMM-Newton obser- 
vations available for the host galaxies of two well-studied 
ULXs (Holmberg II and Holmberg IX), these observa- 
tions allow for the added analysis of Reflection Grat- 
ing Spectrometer (RGS) spectra in addition to spectra 
from the European Photon Imaging Cameras (EPIC). 
The spectral resolution of the EPIC and RGS allow us to 
test different soft component models for the ULX sources, 
to determine the effect of the model on absorption and 
abundance measurements. 

2. SOURCE SELECTION AND DATA REDUCTION 

In iPaper ll we conducted an archival XMM-Newton 
study of ULXs in 32 nearby (< 8 Mpc) galaxies. As 
described in that paper, we extracted spectra for the 
brightest sources in the observations, corresponding to 
> 400 counts. In this study, we chose to further ana- 
lyze the spectra of the objects with the highest num- 
ber of counts (>_5000 counts^). In addition to the 11 
sources from lPaner IL we include an analysis of 3 addi- 
tional sources: the two ULXs in the spiral galaxy NGC 
4559 (observation 0152170501) and the bright source in 
M33 (observation 0102640101). A fuU list of the 14 ULX 



sources, with details of the observations (including expo- 
sure times and count rates), is found in Table 

The two UL X sources in N GC 4559 were origi- 
nally studied W^^gle^^ietsc^^ Bertoldi (1997) using 
ROSAT and 'Crop per et al\ l|2004D with XMM-Newton. 
We follow the naming convention established in these 
papers. Bo th of the se sources (X7 and XIO) were not 
included in iPaper ll because the host galaxy's distance 
is greater than the 8 Mpc limit we initially required. 
However, we include these sources now due to the high 
number of counts (> 5000 counts) in their spectra. Ini- 
tially we did not include M33 X-8 in our ULX sur- 
vey due to its location in the center of its host galaxy. 
This source, however, shows no evidence of being a low- 
luminosity AGN and is more likely a black hole X-ray 
binary (Takano et al. 199^!^). 

Since the initial study of IPaper ll longer exposure 
time XMM-Newton observations have become available 
for three of our sources from the original study. With 
the permission of Tod Strohmayer, we include data from 
the 100 ks XMM-Newton observation (0200980101) of 
Holmberg IX. These data include pn and MOS spec- 
tra of Holmberg IX XMMl and MOS spectra of M81 
XMMl. We also include a 100 ks XMM-Newton obser- 
vation (0200470101) of Holmberg II XMMl that became 
public after the Paper 1 study. The EPIC and RGS spec- 
tra from this ob servation of Holmb erg II XMMl were 
first analyzed bv iGoad et al\ l|2005j) . For Holmberg IX 
XMMl and Holmberg II XMMl, the 100 ks exposures 
provided us with the opportunity to extract and analyze, 
in addition to the EPIC spectra, spectra from the RGS 
detectors. Thus, we include an analysis of RGS spectra 
for both Holmberg II XMMl and Holmberg IX XMMl. 

For the EPIC spectra we added to our original sample, 
we followed the same reduction method as in IPaperll . 
For observations that were processed with an earlier ver- 
sion of the XMM-Newton Science Analysis System (SAS) 
(we used SAS version 6.0^), the observation data files 
(ODE) were used to produce calibrated photon event files 
for the EPIC-MOS and pn cameras using the commands 
emcliain and epchain. The events tables were filtered 
using the standard criteria outlined in the XMM ABC 
Guide. For the MOS data (both MOSl and M0S2 cam- 
eras), good events constitute those with a pulse height 
in the range of 0.2 to 12 keV and event patterns that are 
characterized as 0-12 (single, double, triple, and quadru- 
ple pixel events). For the pn camera, only patterns of 
0-4 (single and double pixel events) are kept, with the 
energy range for the pulse height set between 0.2 and 15 
keV. The selection expression "FLAG == 0" was used 
to exclude bad pixels and events too close to the edges of 
the CCD chips. Time filtering was applied as needed by 
editing the light curve produced in xmmselect. For the 
EPIC observations, time periods in the observation with 
high count rates (flares) in the MOS and pn were cut 
using the command tabgtigen with the 'RATE<' com- 
mand set to 5ctss^^ for MOS detectors and 20ctss~"'^ 
for the pn detector. 

^ In the appendix we show results of simulations to determine the 
number of counts necessary to detect the oxygen and iron absorp- 
tion edges. We found that > 40000 counts are needed to constrain 
iron without having large errors in the measurement. For oxygen, 
> 5000 counts are needed to constrain oxygen without large errors. 
See the appendix for further details. 



High Signal-to-Noise ULX Study 



3 



TABLE 1 

Details of the ULX Sources Analyzed 



Source 


RA (h m 8)=^ 


Dec (o / n)^ 




Obs ID 


Exposure Time (s)"^ 


Count Rate (ctss ^)'^ 


NGC247 XMMl 


nn 
uu 


A7 


uo.o 


-20 47 46.2 


1.54 


0110990301 


3458, 1389, 1379 


0.20, 0.06, 0.06 


NGC253 XMM2 


nn 
uu 


A7 


99 A 


-25 20 55.2 


1.40 


0152020101 


-, 10347, 10304 


-, 0.08, 0.09 


NGC300 XMMl 


nn 
uu 




HQ Q 


-37 42 13.9 


3.11 


0112800101 


6778, 2248, 2453 


0.19, 0.05, 0.05 


M33 X-8 


ni 

UJ- 


OO 


Q 


-1-30 39 36.1 


5.58 


0102640101 


11919, 12147, 12142 


0.06, 0.02, 0.02 


NGC1313 XMM3 


n'^ 
uo 


1 s 

10 


99 ^ 


-66 36 06.2 


4.0 


0106860101 


6960, 2179, 1793 


0.233, 0.08, 0.07 


Holm II XMMl 


08 


19 


28.8 


-1-70 42 20.3 


3.42 


0112520701 


31052, 1257, 10807 


2.72, 0.76, 0.73 














0200470101 


56987,65758, 65766 


3.29, 0.95, 0.94 














0200470101 


41816, 41802 (RGB) 


0.09, 0.11 (RGS) 


M81 XMMl 


09 


55 


32.9 


-1-69 00 34.8 


4.12 


0111800101 


50788, -, 18988 


0.51, -, 0.22 














0200980101 


-, 111910, 111930 


-, 0.111, 0.104 


Hoi IX XMMl 


09 


57 


53.3 


4-69 03 48.7 


4.0 


0112521001 


14976, 6546, 6586 


2.07, 0.64, 0.65 














0200980101 


104010, 111760, 111830 


1.69, 0.51, 0.51 














0200980101 


103840, 103830 (RGS) 


0.03, 0.04 (RGS) 


NGC 4559 X7 


12 


35 


51.8 


-1-27 56 04 


1.51 


0152170501 


34517, -, - 


0.318, -, - 


NGC 4559 XIO 


12 


35 


58.6 


-1-27 57 40.8 


1.51 


0152170501 


34509, -, - 


0.238, -, - 


NGC4631 XMMl 


12 


41 


55.8 


-1-32 32 14 


1.28 


0110900201 


5093, 1969, 1762 


0.13, 0.04, 0.04 


NGC5204 XMMl 


13 


29 


38.5 


-f 58 25 03.6 


1.42 


0142770101 


9981, 3352, 3384 


0.628, 0.177, 0.179 














0142770301 


9231, 2284, 2349 


0.855, 0.247, 0.258 


M83 XMMl 


13 


37 


19.8 


-29 53 49.8 


3.94 


0110910201 


3074, 927, 987 


0.12, 0.033, 0.025 


NGC5408 XMM1<1 


14 


03 


19.8 


-41 22 59.3 


5.73 


0112290601 


5932, 2036, 2077 


0.128, 0.032, 0.033 



^RA and Dec values quoted are the source positions from the EPIC-pn images. 

''Milky Way hydrogen column density along the line of sight from Dickev & Lock mai] 119901) in units of 10^" cm~^ 
'^Exposure times and count rates are listed for the EPIC pn, MOSl, MOS2, and RGS (1 and 2), where available. Note that 

for sources with multiple observations, details for the additional observations are listed below the first observation. Details of 

the specific RGS (RGSl and RGS2) observations are indicated with (RGS). 
'^A second observation of NGC 5408 XMMl is referred to in Section 4.2. This observation (0302900101) with an exposure 

time of 130335 s is proprietary and an analysis will appear in Strohmayer et al. (in prep). 



We extracted source and background spectra along 
with response and ancillary response matrices using 
the SAS task especget. The source spectra were ex- 
tracted from circular regions, typically with radii of 
20 arcseconds. This region was adjusted depending on 
the size of the source and the proximity of the source to 
either another source or the edge of a CCD chip. We ex- 
tracted background spectra from annular regions, except 
when the source was near another source or near the edge 
of a chip. In this case, in order to avoid source confusion, 
we extracted background spectra from a circular region 
located near the source and on the same CCD chip (as 
for Holmberg IX XMMl and M81 XMMl). 

For the RGS spectra we extracted first and sec- 
ond order spectra for the sources Holmberg IX XMMl 
(0200980101) and Holmberg II XMMl (0200470101) us- 
ing the rgsproc command. The RA and Dec values used 
to extract the RGS spectra were obtained from the EPIC 
pn data and are the values quoted in Paper 1. Time fil- 
tering was applied as for the EPIC data with tabgtigen, 
where the 'RATE<' command was set to 0.5ctss^^ for 
Holmberg IX and O.lctss^^ for Holmberg II. Once the 
spectra were obtained, for the RGS as well as EPIC data, 
they were rebinned to require at least 20 counts per bin, 
using the command grppha in LHEASOFT. 

3. SPECTRAL FITTING 

^ Since the original processing of the data in this study, a new 
version of SAS (6.5) had become available. This version includes 
updates to the low energy response files. In the appendix, we show 
results of fitting the pn observations processed with SAS 6.5 (for 
pn counts > 5000). We found that the model parameters are not 
significantly different between the SAS 6.0 EPIC spectra and SAS 
6.5 pn spectra. 



Spectral fitting proceeded using XSPEC vl 1.3.1. For 
the RGS spectra, we simultaneously fit the first order 
spectra from both RGSl and RGS2 in the RGS band 
(0.33 - 2.5 keV). For the EPIC spectra, we fit the pn 
and MOS spectra simultaneously in the 0.3-10 keV en- 
ergy range. We allowed a free normalization constant to 
account for the differences in flux calibration between the 
three EPIC cameras. Both the RGS and EPIC spectra 
were fit s eparately . 

In the jPaper ll study, we had fit all of the sources 
with three standard models: an absorbed power law, 
an absorbed bremsstrahlung model, and an absorbed 
combined blackbody and power law model. We used 
the XSPEC model wabs to account for absorption from 
the Milky Way and the host galaxy/ULX contribution. 
This model is a ph oto-electric a bsorption model using 
the cross-sections of lMorri^ 9n ^ McCammon ( 1983) and 
the solar abundances of lAnders fc Ebiharal <|1982|) . We 
found that the spectra of the brightest ULXs were typ- 
ically best-fit by an absorbed combined blackbody and 
power law model (fit in XSPEC as wabs*wabs* (bbody 
+ pow) where the first wabs model was fit to the 
iDickev fc LockmanI l)1990() Milky Way value and the 
second wabs model was fit to the remaining host 
galaxy/ULX contribution). Likewise, we began our 
study by fitting the additional sources (M33 X-8, NGC 
4559 X-7, and NGC 4559 X-10) with the same three mod- 
els noted above. We found that these sources were well 
fit by the absorbed blackbody and power law model with 
xVdof « 1.0. 

For M33 X-8, the source was well-fit with an absorbed 
blackbody and power law with best-fit parameters: uh = 
l.Qlt'.Z X 10^^ cm-2, kT= Q.lAt','^ keV, P = 2.46^;;:°^, 
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TABLE 2 

Spectral Fits for EPIC spectra with tbabs*tbvarabs*edge*(grad + pow) model 



Source 






abund.'^ 


Fc abund.*^ 


iviass l^i-Vi0 J 


A/f 1 aV c 




1 




counts'^ 


Holmll XMMl 


0, 


-,2+0-01 
— 0.01 


0, 


■^'J-o.ii 


0, 


n+0.15 
■'J-0.0 


141+21 

— oz 


0.15+^9 
—U.uy 


2, 


og+0.05 
—0.05 


1.11 


342874 




0, 


7+0. 02 


1, 


07+0.23 

-0.42 


3, 


15+1.29 


791+209 


0.07IH3 


2, 


29 + g°Q 


0.98 


43116 


HolmIX XMMl 


0, 


]^g+0.02 


1, 


O/I+0.12 

■'-'^-o.is 


2, 


07+0.12 
■'^'-0.13 


382t?6 


O.OQlH' 


1, 


3g+.025 


1.05 


148061 




0, 


oQ+0.08 
— u.uy 


1, 


Qo+0.29 
■'-'°-0.40 


3, 


ni+1-09 
■^^-1.86 


1181+1^16 


— u.u / 


1, 


7g + 0.09 

— 0.09 


0.64 


28108 


M33 X-8 


0, 


2Q+0.04 


1, 


(,.+0.13 
■'-'*-0.18 


1, 


74+0-62 
■'^-0.74 




1 n2+''-°'' 

-L ■ Qg 


2, 


g2+0.34 


0.97 


123903 


M81 XMMl 


0, 




1, 


on+0-08 

■^"-0.10 


2, 


04+0-44 
■'J*-0.57 


8.6ltfof 


9 09+44.8 

. 0.^ -j^ g-j^ 


4, 




1.02 


69776 




0, 


27+0.21 


1, 


49+0.58 
■^^-0.50 


1, 


g7+1.61 
■°'-1.67 




01 +4.36 

O.Ol_g gg 


2, 




0.91 


31731 


NGC253 XMM2 


0, 




1, 


c:c+0.20 
■'''^-0.29 


2, 


70+1.17 
■ '•-'-1.46 




0-05+^3 

— U.Uo 


2, 


22+0.10 


0.96 


20651 


NGC5204 XMMl 


0, 


■10^003 


1, 


49+0.55 
■^^-0.75 


0. 


0+1.84 
■"-0.0 


4641366 


0-041^2 


1, 


92+ggg 


0.96 


16717 




0, 


-|^2+0.05 
— 0.03 


0, 


77+0.60 
■ ' ' -0.77 


0, 


0+1.70 


449t^90 


— U.Uo 


2, 


■O^-o.ll 


0.93 


13864 


NGC1313 XMM3 


0, 




1, 


07+0.13 
"J' -0.14 


0, 


0+0.34 
■^-0.0 




0-08^0.02 


2, 


■66IH9 


1.02 


10932 


NGC300 XMMl 


0, 




2, 


■^"-0.43 


0, 


0+1.77 
■^-0.0 


4171^75 


0.02l«:« 


2, 


.46i«:i; 


1.01 


11479 


IN 4ooy A- ( 


U. 


-,0+0.05 
■J^"-0.03 


0, 


■^°-0-46 


0, 


0+1.00 
■"-0.0 


7cc:+901 
l'JJ_443 


nfi+0.12 


z. 


-10+0.10 


U.oo 


iZOUD 


NGC4631 XMMl 


0, 


90+0. 06 
■^"-0.06 


0, 


C9+0-43 
■"^-0-58 


0, 


1 1 +1.86 
■^^-0.11 


c c: + 11.4 
''•'^-IS.O 


0-841H2 


5, 


7-1 +2.03 
■'^-0.31 


1.06 


8824 


NGC5408 X-1 


0, 


„Q+0.03 
■'-'^-0.03 


1, 


Qg + 0.26 

■^='-0.52 


5, 


0+0.0 
■'J-1.68 


14771^7? 


"■l"-0.09 


2, 


c- .1+0. 20 
■•''^-0.16 


0.96 


10045 


NGC4559 X-10 


0, 


-, 0+0.05 
■-'-^-0.04 


1, 


O0+O.6I 
■^°-0.81 


0, 


0+2. 30 
■"-0.0 


5 04+7.88 


2 74+2.52 
^- '^-1.86 


2, 


09+0.52 
■"='-0.28 


0.77 


8837 


NGC247 XMMl 


0, 


44+0-18 
■^^-0.12 


1, 


■'-'='-0.28 


0, 


0+1.45 
■"-0.0 


1 71 7+3205 
-l'-L'_929 


n 05+0.02 
'^■'-'0-0.03 


3, 


gO+6.40 
■""-6.60 


0.60 


6226 


M83 XMMl 


0, 


1 0+0. 09 
■^'-'-0.10 


1, 


7„+1.12 
■"^-1.65 


0, 


0+3.92 
■"-0.0 


6.341133* 


Q2+0-™ 
"•'^^-0.70 


2, 


fi4+0.88 
■"^-0.64 


0.83 


4988 



^ Hydrogen column density d etermined from tbvarabs in units of 10^^ cm ^. The Galactic value of n^f was fixed to 
the lDickev fc LockmanI 1199(1) value with the tbabs model. 
''Element abundance relative to the Wilms solar abundance from the tbvarabs model 

Ratio of mass accretion rate from the grad model to Eddington accretion rate (see Section 4) 
'^Total number of photon counts from pn and MOS detectors 



TABLE 3 

Spectral Fits for RGS spectra with tbabs*tbvarabs*edge*(grad + pow) model 



Source 




abund.'' 


Fc abund.*" 


Mass (M0) 


M/MEdd" 


F 


xVdof 


counts'* 


Holmll XMMl 


0.08«;«^ 


0-65lE;;«i 


00+1.51 


1 47+8® 
-■-^'-83 


24+°-" 
"•^^-0.12 


1-58113 


444.4/442 


9521 


HolmIX XMMl 


0-29trol 


0-68l»;30 


.ig+1.36 
"•■^"-0.36 


774+4226 
' '^-570 


"■"^-0.03 


1 43+0.45 

-'-•^•-'-0.98 


290.4/339 


10807 



^Hydrogen column density determined from tbvarabs in units of 10^^ cm ^. The Galactic value of xih was fixed to 
the Dickcv & Lockman 1 1990) value with the tbabs model. 
''Element abundance relative to the Wilms solar abundance 

Ratio of mass accretion rate from the grad model to Eddington accretion rate (see Section 4) 
''Total number of photon counts from RGSl and R,GS2 detectors 



and = 1579.8/1533 dof, where nn represents the host 
galaxy/ULX hydrogen column (the Milky Way contri- 
bution was fixed to thelDickev fc LockmanI ljl99(]|) value 
listed in Table In lPaoer ll we noted that ULX sources 



well-fit with the combined blackbody and power law 
model with a higher disk temperature (w 1 keV) and 
a lower flux were often well-fit by an absorbed Comp- 
tonization model (XSPEC model compst). We note 
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Fig. 1.— EPIC spectrum of observation 0200980101 of Holm- 
berg IX XMMl. This plot shows the spectrum fit with the 
tbvarabs*tbabs*edge* (grad + pow) model and the residuals from 
this fit. The best fit parameters are listed in Table |51 



that for M33 X-8 the inverse Compton scattering model 
also fits the data well (wabs*compst) but with a = 
1614.9/1536 dof. Despite the larger value, the compst 
model better fits the residuals from the spectra's sloping 
high energy tail. The best-fit absorbed blackbody and 
power law parameters for NGC 4559 X7 (12 h 35 m 51.8 s, 
27° 56' 4") were: hh = l-^ltom x lO^^cm-^, kT= 
0.13;[!-^2keV, r = 2.l6to M, and = 410.3/369 dof. 
For NGC 4559 XIO (12 h 35 m 58.6 s, 27° 57' 40.8"), we 
found: TiH = 1.14to;oi x lO^^cm'^, kT= 0.96l°;26 keV, 
r = 2.131^-^^ and = 292.3/330 dof. XIO, as noted in 
ICroDPer et all lj2004j) . is also well fit by a Comptoniza- 
tion model. 

For the sources M33 X-8 and NGC 4559 XIO, we noted 
that their spectra were well-fit by either an absorbed 
blackbody and power law or an absorbed Comptoniza- 
tion model. This brings up one issue surrounding ULX 
spectra. Namely, ambiguity over their spectral form. In 
this paper, we assume that there is an analogy between 
ULXs and Galactic black hole binaries. Thus, we assume 
that the accurate ULX spectrum for these bright sources 
is a hard component (well modeled by a power law) and 
a soft component (which we assume as a thermal con- 
tribution from an accretion disk). The form of the soft 
component, in particular, will affect the measured ab- 
sorption and abundance values. Thus, an investigation 
of this component and its affect on the absorption model 
is important. 

The soft component of ULXs is most often modeled 
as a thermal component originating from an accretion 
disk surrounding a central black hole. There are nu- 
merous disk models applied to model this possibly ther- 



mal component. In iPaoer ll we modeled this compo- 
nent as a simple blackbody (in XSPEC bbody). While 
an accretion disk is expected to have a range of tem- 
peratures, empirically a single blackbody is a good fit 
to low signal-to-noise spectra. A simple absorbed black- 
body and power law model was used for Galactic black 
hole X-ray binaries in the 1980s when the quality of 
data for these sources was analogous to that for ULX 
sources today. As the next step in accretion disk models, 
the soft component of ULX spectra is often modeled as 
an optically thick, geometrically thin, multi-component 
blackbody disk (a m ulti-component disk or MCD model, 
diskbb in XSPEC) lIMitsuda p.t ffi.iri98l . 

Two disk models that are used to fit the soft com- 
ponent of ULX spectra with more physical accuracy 
are the XSPEC models diskpn and grad. The diskpn 
model is an extension of the MCD (diskbb) model which 
uses a ps eudo-Newtonian p otential. The grad m odel 
l)Hanawal [T'989: Ebisawa. Mit suda. fc Hanawalll99lD is a 
multi-component disk model that unlike the diskbb or 
diskpn model, incorporates the affects of general rela- 
tivity. One of the advantages of the grad model is that 
it fits the spectra for mass (Mgrad) and mass accretion 
rates (M) given a few initial assumptions (distance to 
the source, disk inclination angle, and the ratio of the 
color temperature to the effective temperature). 

In addition to thermal models, other models have 
been suggested to explain the soft component. Two 
of these models are ionized refiection and the warm 
absorber model. Both of these models have been ap- 
plied to low-redshift PC Quasars, sources with black- 
body temperatures of 150 eV. This is relevant to ULXs 
since many spectral fits of ULXs require cool acc retion 
disk temperatures of approximately 100 eV (Millcr ~^ all 
'2003; Miller. Fabian & Miller 2004; Roberts et al. 2005; 
■Winter. Mushotzkv. fc Revno lds 2005). The reflection 
model suggests that the soft component results from X- 
ray ionized reflection. In this model, back-scattering 
and fluorescence of X-rays in the disk, as well as radia- 
tive recombination, cause elements with smaller ioniza- 
tion potentials (e.g. C, O, N) to become highly ionized. 
iR.oss fc FabianI (2005) note that a relativistically blurred 
X-ray ionization model folded through an XMM-Newton 
pn response matrix, is well-fit by a blackbody with a 
temperature of 150 eV, the same value that is seen in PC 
Quasars and many ULXs. In the warm absorber model, 
absorption edges and lines from an absorbing material 
close to the X-ray source appear as a thermal compo- 
nent in poorer quality spectra. With increased spectral 
resolution, the numerous absorption edges and lines are 
distinguishable. The warm absorber may be the result 
of a stron g, mildly relativistic wind from the disk as sug- 
gested bv IGierliiiski fc Done! l|2004D . Indeed, the pres- 
ence of a warm absorber is well known and studied in 
many Seyfert galaxies (e.g. NGC 3783 and MCG-6-30- 
15). 

As noted above, for subsequent spectral fits we as- 
sume a thermal model for the soft component of the 
ULX spectra. In the following section (Section 4), we 
will discuss the effect different thermal models have 
on the measured absorption values. In order to mea- 
sure the hydrogen column density, we fit the spec- 
tra of sources listed in Table with the more so- 
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phisticated T uebingen-Boulder IS M absorption model 
of Wilms. All en, fc McCravl l|200(l (tbabs, tbvarabs in 
XSPEC). This model accounts for X-ray absorption re- 
sulting from contributions from X-ray absorption from 
the gas phase of the ISM, grains, and molecules. The 
model uses updated solar abundances and photoioniza- 
tion cross-sections. For both the EPIC and RGS spec- 
tra, we accounted for Galactic hydrogen absorption by 
setting the column density of the tbabs model equal to 
the Milky Way hydrogen column density along the line 
of sight to the host galaxy . These column densities are 
quoted in Table 1 and are all less than 6 x 10^°cm~^. 

Since the measured X-ray column densities for ULXs 
are typically an order o f magnitude higher than the 
Dickey fc LockmanI ljl990(l Milky Way values (Fig. 8 of 
Paperlfi . we are confident that the additional absorp- 
tion measured is not from the Milky Way. To deter- 
mine the host galaxy's hydrogen column density and the 
abundances of elements along the line of sight (oxygen 
and iron), we used the tbvarabs model. The tbvarabs 
model accounts for X-ray absorption due to photo- 
ionization. It includes the effects due to the II2 molecule 
and depletion of metals in grains. The model allows for 
individual fits to abundant elements (He through Ni) , H2 , 
and depletion of elements in grains. We initially allowed 
the hydrogen column density to vary, fixing all other 
tbvarabs parameters to their defaults. The Galactic 
column density remained fixed (using the tbabs model) , 
while the other parameters (power law and blackbody 
components) were allowed to vary. This was modeled 
as tbabs*tbvarabs* (bbody + pow) in XSPEC, where 
tbabs was fixed to the Milky Way value and tbvarabs 
was used to fit the absorption from the ULX/host galaxy 
contribution. With the best-fit hydrogen column density 
(from tbvarabs), we allowed the oxygen abundance and 
then the iron abundance to vary from the solar abun- 
dances. Allowing these parameters to float provides mea- 
surements of the depth of the oxygen K-shell edge at 
542 eV a nd the iron L-shell edge at 8 5 1 eV. 

As in IBaumgartner fc Mushotzkvl lj2006D . we found 
that for our 14 sources, the oxygen absorption values 
from the tbvarabs model yielded different values for the 
EPIC pn and MPS spectra. Thus, we foll ow the proce- 
dure of IBaumgartner fc Mushotzkvl ||200^ in adding an 
edge model to account for the differences. We add an ex- 
tra edge component to the MOSl and M0S2 detectors at 
an energy of 0.53 keV with optical depths of 0.22 and 0.20 
respectively (modeled as tbabs*tbvarabs*edge* (bbody 
+ pow))'^. This is the template model used for all spec- 
tral fits mentioned throughout the rest of the paper, 
where other thermal models are used in place of the sim- 
ple blackbody (bbody) model as indicated. 

An additional problem we note involves the hy- 
drogen column densities along the line of sight. 
IBaumgartner fc Mushotzlcvl l)2006|) note that a column 
density (uh) greater than 0.5 x 10^"'^ cm~^ is necessary 
to allow for good oxygen abundance measurements. We 
also note that at very high hh values the oxygen edge 
is undetectable, where high absorption values cause the 
signal-to-noise to be too low to detect the edge. 

^ T he edge parameters used by Baumgar tner Hi Mushotzkvl 
120061) were calculated for differences in calibration using SAS ver- 
sion 5.4. In the appendix, we show that the pn values using SAS 



4. NATURE OF THE SOFT COMPONENT 

In the previous section we outlined three possible ori- 
gins for the soft component in ULXs (X-ray reflection, 
thermal emission from an accretion disk, and a warm ab- 
sorber model). For our purposes, we assume a thermal 
component origin. We chose the thermal model since, as 
described in Paper I, the spectra of ULXs are well de- 
scribed with models used to fit Galactic black hole X-ray 
binaries. To gauge the effect different disk models (i.e. 
bbody, diskbb, diskpn, grad) may have on the measured 
abundances, we discuss different disk models in section 
4.1. In section 4.3, we assess the usefulness of the model 
(absorbed disk and power law) for measuring the galac- 
tic/ULX hydrogen column density and the oxygen and 
iron abundances by discussing the physical plausibility 
of this model. 

For sources with at least 5000 counts, the 14 
sources listed in Table ^ we fitted the grad model 
for the soft-component of the spectra (XSPEC model 
tbabs*tbvarabs*edge* (grad + pow)). We assumed 
a disk inclination angle of 60 degrees and a ratio 
Tcoior/Te// = 1.7 (the default value, typical of stellar 
mass BH sources) . In Table |21 we list the best-fit pa- 
rameters for EPIC spectra using this model. Rows with 
no source name indicated represent an additional obser- 
vation of the previous ULX (as indicated in Table 
Table O provides the best-fit parameters for the RGS 
spectra. The important measurements to note are the 
host galaxy/ULX contribution to the hydrogen column 
density, oxygen abundance, and iron abundance from the 
tbvarabs model. All errors quoted in this paper are 
for the 90% confidence level for one degree of freedom 
(Ax^ — 2.76). We provide a representative spectral fit in 
Figure^for the long observation of Holmberg IX XMMl. 

4.1. Comparison of Thermal Disk Models 

We chose the general relativistic disk model, grad, be- 
cause it is the most physically accurate of the various sim- 
ple accretion disk models. Also, the grad model requires 
few initial parameters while making a direct calculation 
of the mass and mass accretion rate of the black hole. 
In order to compare the grad results with an alternate 
model, we fit the highest signal-to-noise observation of 
Holmberg IX XMMl with the diskpn model (with the 
inner radius of the disk set to the radius of marginal 
stability or 6 times the Schwarzschild radius) in place 
of the grad model. We found that the values obtained 
agreed with those of the grad model. The mass from 
the diskpn model was slightly lower at 358 M© (cor- 
responding to kT= 0.23 keV) and the iron abundance 
was slightly higher at 2.33 compared to the grad model. 
Since the diskpn and grad models yield similar values, 
either of these disk models would be sufficient for de- 
termining the hydrogen column density and oxygen and 
iron abundances. 

In order to test the effect the thermal model has 
on the tbvarabs measured parameters (i.e. hh and 
abundances), we wanted to compare the results from 
the more physical disk model (grad) with less phys- 
ical models (bbody and diskbb). In order to com- 

6.5 are consistent with SAS 6.0 pn-|-MOS derived values. Thus, 
the correction quoted is valid for our purposes with SAS version 
6.0. 
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TABLE 4 

Alternate Thermal Model Spectral Fits to ULX Sources with > 20000 Counts 



Source ii//'' Oxygen abundance'' Iron abundance'' kT T x^/dof counts'^ 



tbabs*tbvarabs*edge* (bbody + pow) 



Holmll XMMl 


0. 


— 0.01 


0, 


.97+°-^^ 


0, 


Q+0.32 


0. 


24+0.01 
— 0.01 


2, 


4]^+0.03 
— 0.03 


2216.5/2018 


342874 




0. 


14+0.03 
-0.02 


1, 


07+0.38 

-0.39 


4, 


00+0. 62 
■■^^-1.8 


0. 


lg+0.01 
^"-0.01 


2, 


01 +0.08 
"^-^-0.04 


952.7/973 


43116 


HolmIX XMMl 


0. 


14+0.02 
-^^-0.01 


1, 


CI +0.19 
■3-^-0.15 


3, 


07+0.99 
"^'-0.80 


0. 


1Q + 0.01 

^^-0.01 


1, 


40+0.02 

■^'J-o.oi 


2856.7/2752 


148061 




0. 


r,Q + 0.10 

^•5-0.04 


1, 


07+0.43 
■o'_0.44 


4, 


7C+0.25 
■ '"^-1.94 


0. 


1 7+0.02 
-^'-0.03 


1, 


71 +0.14 
■"-0.06 


559.5/880 


28108 


M33 X-8 


0. 


99+0.01 
^^-0.14 


1, 


94+0.10 
■^^-0.10 


2, 


00+0.31 
"^"-0.60 


0. 


7(5+0.02 
'"-0.02 


2, 


CI +0.07 
■•^-^-0.06 


1528.8/1533 


123903 


M81 XMMl 


0. 


oQ+0.03 
■^^-0.02 


1, 


1 O+0.15 
■-'-•^-0.07 


1, 


70+0.65 
■ '°-0.40 


0. 


gQ + 0.03 

="-'-0.02 


2, 


7n+0.05 
■ "^-0.05 


1240.2/1241 


69776 




0. 


9fi+0.09 
^°-0.04 


1, 


46+0.26 
■^"-0.49 


1, 


fii +2.02 
■°'-1.61 


0. 


99+0.38 
="^-0.28 


1, 


7f;+0.13 
■ '"-0.09 


627.4/668 


31731 


NGC253 XMM2 


0. 


99+0.05 
■^^-0.05 


1, 


91 +0.40 
■^-■--0.37 


1, 


cq+1.67 
■■^•^-1.53 


0. 


70+0.14 
'■""-O.ll 


2, 


1C+O.O8 
■^^"-0.23 


460.8/496 


20651 



tbabs*tbvarabs*edge* (dlskbb + pow) 



Holmll XMMl 


0, 


, , +0.001 

■'1-0.01 


0, 


fic+0.09 
■°3-0.05 


0.0+0.06 


0, 


01+0.02 
■'-'^-0.02 


9 00+0. 05 
^■•-"'-0.04 


2241.8/2018 


342874 




0, 


1 c+0.03 
■'■^-0.02 


1, 


nc+0.34 
■^^^-0.38 


2.67ti;?i 


0, 


91 +0.04 
■^1-0.01 


9 97+0.06 
-0.07 


954.2/973 


43116 


HolmIX XMMl 


0, 


1Q+O.O2 
■1^-0.02 


1, 


00+0.11 
■■^"-0.16 


9 oQ+0.58 
^•■^^-0.84 


0, 


91+0.02 
■^^-0.01 


1 09+0.02 


2878.9/2752 


148061 




0, 


OQ + 0.11 

■'-'"-0.07 


1, 


of;+0.32 
■'-'"-0.38 




0, 


9n+0.04 
■^"-0.03 


1 79+0.10 
1- '^-0.07 


559.8/880 


28108 


M33 X-8 


0, 


1Q+0.04 

■i-^-o.os 


1, 


04+0.15 
■"^-0.16 


1 7C+O.68 
l- '■'-0.71 


1, 


1 7+0.05 
■"-0.05 


2.47«:?i 


1530.4/1533 


123903 


M81 XMMl 


0, 


Rq+0.08 

■"•^-o.is 


1, 


1Q+0.09 
■1^-0.15 


^•"'-0.69 


1, 


00+0. 03 
■■^"-0.08 


4.08«;r9 


1247.7/1241 


69776 




0, 


97+0.18 
■^'-0.08 


1, 


42+0.57 
■^^-0.49 


1 63+'-'^'* 
^•°'^-1.63 


2, 


Qg+0.55 
■""-0.97 


1 98+'-'"' 
1-^^-0.58 


626.9/668 


31731 


NGC253 XMM2 


0, 


91 +0.10 
■^'-0.07 


1, 


10+0.53 
■'"-0.37 


1 64+'-'^* 


1, 


97+0.22 
■^'-0.27 


^•^'-0.62 


461.3/496 


20651 



"Hydrogen column densi ty determined from tbvarabs in units of lO'^'^ cm ^. The Galactic value of njj was fixed to the 
IDickev LockmanI 1199(1^ value with the tbabs model. 
''Element abundance relative to the Wilms solar abundance from the tbvarabs model 
Total number of photon counts from pn and MOS detectors 



pare the results from these three models, we chose 
to examine the best-fit parameters for spectra with 
the highest signal-to-noise (using sources with at least 
20000 counts). Thus, in Table H we hst the best- 
fit parameters obtained using an absorbed black- 
body and power law (tbabs*tbvarabs*edge* (bbody + 
pow)) and an absorbed diskbb and power law model 
(tbabs*tbvarabs*edge* (diskbb + pow)) for sources 
with at least 20000 counts. Subsequently in this paper, 
when we refer to the bbody, diskbb, or grad models, we 
are referring to the model fits used in Tables El 

The curvature, or shape of the spectra at low ener- 
gies, for the bbody, diskbb, and grad models is different, 
presumably affecting the column density and absorption 
values. From a comparison of these models, compar- 
ing the mean values from all of the observations with 
> 20000 counts hsted in Tables El and H we find that 
the column densities obtained with the diskbb model 
are nearly identical to those with the grad model (the 
exception is the source NGC 253 XMM2). The MCD 



(diskbb) model hydrogen column densities are approx- 
imately 6% higher than those of the blackbody model 
(excluding the first M81 XMMl observation). The aver- 
age power law indices for the three models agree within 
a factor of 5%, with the bbody model having the lowest 
and the grad having the highest average power law index 
(excluding the first M81 XMMl observation). Thus, the 
difference in both power law index and hydrogen column 
density is negligible on average, with the exceptions of 
NGC 253 XMM2 and M81 XMMl. 

For a comparison of the derived abundances, the oxy- 
gen abundance is roughly the same between the grad 
and bbody model. These values are approximately 
10% higher than those derived from the MCD (diskbb) 
model. The iron abundances, however, vary more from 
model to model. The blackbody model average iron 
abundance is w 13% higher than that of the grad model 
and w 18% higher than that of the diskbb model. Thus, 
there are no significant changes to the model parameters 
in changing the base accretion disk model. These results 
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channel energy (keV) 

Fig. 2.— Plot of the EPIC unfolded spectrum for the NGC 
5408 XMMl ULX with the high mass solution (see Table 2). With 
the spectrum, the absorbed grad and absorbed power law model 
components are plotted with light gray lines while the combined 
model is plotted in dark gray. The thermal component clearly 
dominates the low energy spectrum with this model. 




channel energy (keV) 

Fig. 3.— Plot of the EPIC unfolded spectrum for the NGC 
5408 XMMl ULX with the low mass solution (see Table 5). With 
the spectrum, the absorbed grad and absorbed power law model 
components are plotted with light gray lines while the combined 
model is plotted in dark gray. Notice that the low energy spectrum 
is dominated by the power law component with a weak contribution 
from the disk model (grad). 



are verified by simulations described in the appendix. 

4.2. Degenerate Model 

In modeling our sources with the 
tbvarabs*tbabs*edge* (grad+pow) model, we found 
that some of the ULX sources were well-fit by a model 
where the power law comp onent dominates the low en- 
ergy spectrum, also seen in lStobbart. Roberts, fc Wilms' 
pOOfi) . This type of spectral fit typically yields a steeper 
power law index (F > 3.0) and a low mass (M< 10 Mq). 
For two sources, M33 X-8 and M81 XMMl, this model 
fit was a much better fit than a higher mass model 
with Ax^ of 476 and 190 respectively. However, for 
some sources, there was a degeneracy between the two 
models (high mass and low mass/steep power law). We 
illustrate this degeneracy with observation 0112290601 
of the source NGC 5408 XMMl, showing the high mass 
model and spectral fit in Figure [3 and the low mass 
model and spectral fit in FigureOl For spectra exhibiting 
this degenerate solution, we include the low-mass/steep 
power law fits in Table These fits all exhibit, in 
addit ion to low masses, solutions w i th M >> MEdd- 

In iStobbart. Roberts, fc Wilms! (|2006), the authors 
discussed the same issues in fitting the XMM spectra 
of 13 ULXs. They found that two sources, M33 X-8 and 
NGC 2403 X-1, were best fit by the model with a power 
law fit to the low energy portion and thermal model at 
higher energy. They also indicated six sources where 
an ambiguity existed between the two models. To un- 
derstand the spectra of the sources where both models 



provided good fits to the data, we further investigated 
the spectra of sources with multiple observations. 

For Holmberg II X-1 and Holmberg IX X-1 our spec- 
tral fits include multiple observations (a shorter and a 
100 kilo second observation). For each of these sources, 
we found that the shorter observation could be fit with 
either a high mass or low mass solution. When we fit 
the 100 kilo second observation, however, the high mass 
model was a much better fit. The Ax^ values between 
the high mass and low mass solutions for the kilo second 
observations were 135 and 88, respectively. To test this 
further, we also fit a 100 kilo second observation of NGC 
5408 XMMl with both of these models. This observation 
(0302900101) is a proprietary observation whose spectral 
and temporal analysis will appear in Strohmayer et at 
(in prep). We processed the pn data with SAS 6.5, follow- 
ing the same procedure as noted in the Data Reduction 
section. Fitting the spectrum with both models (high 
mass and low mass solution) we found a Ax^ value of 
170, favoring the high mass model. 

Thus, we find that for sources that are well fit by either 
model (showing a degenerate solution of either high mass 
or low mass/steep power law) the high mass solution is 
the best fit when a higher count spectrum is obtained 
(as for Holmberg IX XMMl, Holmberg II XMMl, and 
NGC 5408 X-1). Though we fist the alternate model 
column density and absorption values in Table 13 we use 
the parameters listed in Table |21 throughout the paper 
(the high mass solutions). 
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TABLE 5 

Spectral Fits for EPIC spectra with tbabs*tbvarabs*edge*(grad + pow) model for Sources with 
Degenerate Solution (Low Mass/Steep Power Law) 



Source 






abund.'' 


Fe abund.*^ 


Mass (M0) 


M/MEdd" 


r 


X 


2/dof 


counts'* 


Holmll XMMl 





nn4-0 04 
2OI0.M 


n OQ+0 17 


+0.83 
^•^^-1.11 


9 


94+1.52 
^^-1.24 




1 n + 31 





99 


43116 


TT 1 TAT' AT" TV TTV T-1 

HolmIX XMMl 





OQ+0.13 

■^^-0.12 


1 r,c+0.21 
^^■'-'0-0.33 


r, no + 0.95 
2.08+145 


3 


01+1.47 
y^-1.59 


99 0+1.79 
^Z.U_2 gg 


„ +0.79 

■^■0^^-0.87 





64 


28108 


NGC253 XMM2 





10+0. 17 
-■-" — O.OS 


1 10+°-'^' 


1 or, + 1.90 

— 1.32 


6 


09+4.19 
°^-2.22 


1 97+0.71 
—0.40 


9 1 c+0.84 
^- .'-'-'_0.58 





93 


20651 


NGC5204 XMMl 





10+0. 05 
-'-°-0.04 


1 nn+°'3° 


^■'^-0.0 


3 


91 +0.98 
^-^--0.95 


fi oR+0.67 
D.»b_o.09 


q 90+0.39 

■J-^°-0.38 





95 


16717 







9n+006 
^'-'-0.05 


fi'5+0 35 
U.DJ_Q 4g 


^•'J-0.0 


2 


05+2.13 


Y 07+1.35 
' -2.08 







94 


13864 


NGC300 XMMl 





r,i +0.05 
^^-0.03 


1 5Q+O.22 


^■'^-0.0 


1 


fjq+0.85 
0^-0.32 


n 70+0.18 

'J- '■^-0.13 


3.89t»:^^ 


1 


00 


11479 


N4559 X-7 





9n+0.06 
^'^-O.OS 


0.341":^ 


^•'-'-0.0 


5 


'7Q+2.82 
'^-2.75 


c CQ + 1.44 
6.63_2 42 


q 1 n + 0. 55 

'-'■^'^-0.50 





83 


12506 


NGC5408 X-1 





10+0. 05 
-^-"-0.04 


0.87tEJ:2 


Q 45+1.75 
'-'■^^-0.45 


2 


90+1.57 
^^-1.28 


^■'''-1.02 


^■'-'^-0.34 





89 


10045 



^Hydrogen column density determined from tbvarabs in units of 10^^ cm ^. The Galactic value of n^f was fixed to 
the Dickey & Lockman 1 1990) value with the tbabs model. 
''Element abundance relative to the Wilms solar abundance from the tbvarabs model 

Ratio of mass accretion rate from the grad model to Eddington accretion rate (see Section 4) 
'^Total number of photon counts from pn and MOS detectors 



As noted, M33 X-8 and M81 XMMl, sources with very 
high number of counts, were not well-fit with a standard 
disk at low energy, power law at high energy model. They 
were best fit with the steep low energy power law and hot 
disk model shown in Figure |3| Along with these sources, 
NGC 4559 X-10 and M83 XMMl were also well-fit by 
this model. We will further discuss these sources in the 
following subsection. 

4.3. Physical Plausibility of the Accretion Disk models 

In Paper 1 , we had found that ULX spectra are consis- 
tent with the high/soft and low/hard states of Galactic 
black holes. We had classified the sources in this study 
as consistent with the high/soft state. The high/soft 
state, as stated earlier, is characterized by emission from 
an accretion disk and a Comptonized power law tail. 
In order to investigate whether the mass and accre- 
tion rate results from the grad model make physical 
sense in terms of the X-ray binary model, we present 
a comparison of the derived accretion rate in Edding- 
ton accretion rate versus the mass in Figure 01 The 
parameter MEdd was computed as Msdd = 

TjC 

where ry, the efficiency factor, was set to 0.06 and M 
is the value M-grad- We find that for the sources with 
Mgrad > 100 Mq, the accretion rate is computed to 
be below 40% of the Eddington rate. This is assum- 
ing, as the grad model assumes, a Schwarzschild black 
hole. Noting that M/MEdd is equivalent to L/L^dd 
for most disk solutions, the L / l^Edd values for the 
sources with Mg^ad > 100 M© are consistent with those 
of Galactic black hole X-ray binaries in the high state 
(L/L Rrfrf 0.05-1.0) ijMcClintock. Naravan. &: Rvbickil 
I2004D . Thus, the sources with Mg^ad > 100 Mg do have 
accretion rates that are predicted from scaling up (in 
mass) observed high state Galactic black holes. In addi- 
tion, the spectral fit parameters for NGC 4631 XMMl, 
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Fig. 4. — Black hole mass versus mass accretion rate (per Ed- 
dington accretion rate) obtained from the general relativistic disk 
model (grad). The sources with model parameters indicating high 
masses correspond to accretion rates that are below 40% of the Ed- 
dington accretion rate for the given mass. In Paper 1, we noted a 
temperature gap in the range of 0.26 keV to 0.50 keV from spectral 
fits utilizing an absorbed blackbody and power law model. In this 
figure, we see that there is a gap in mass from the parameters of 
the absorbed general relativistic disk and power law model. The 
values shown are results for the EPIC spectra only, as recorded in 
Table 2 (for all sources in Table 2). Therefore, there are multiple 
points for the sources with multiple observations. We plotted these 
sources (Holmberg II XMMl, Holmberg IX XMMl, M81 XMMl, 
and NGC 5204 XMMl) twice to confirm that their masses, as de- 
termined by the grad model, do not vary. 



with an estimated mass of 5.5 Mq and an accretion rate 
L/Lsrffj of 0.84, are also consistent with the standard 
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high state Galactic black hole model. This source is likely 
a normal stellar mass black hole X-ray binary in an ex- 
ternal galaxy. 

The remaining sources with Mgrad << 100 (M33 
X-8, M81 XMMl, NGC 4559 X-10, and M83 XMMl) 
yielded L/LEdd ratios in the range of 1 - 3. These 
sources are also those described in the previous section 
where the power law component fits the low energy spec- 
trum. They are also well fit by a Comptonization model 
and correspo nd to a sub-class of high luminosity ULXs 
described in .Paper li Due to the luminosity and mod- 
eled disk temperature (kT IkeV), we suggested that 
these sources were very high state stellar mass black hole 
systems. In order to be consistent with the black hole ac- 
cretion model assumed in this study, the spectra of these 
sources should be the result of Comptonization from a 
thermal disk spectrum. 

To test this further, we fit the spectrum of the 
highest count source of this type (M33 X-8) with 
an absorbed thermal disk and Comptonization model 
(tbabs*tbabs*(diskpn -f comptt)). The comptt model 
has the parameters: a seed temperature (keV), a plasma 
temperature (keV), and optical depth of the medium. 
We used the diskpn model in place of the grad model 
since the former provides a disk temperature to which the 
comptt seed temperature can be fixed. This provides a 
physical model, where the thermal disk supplies the en- 
ergy for the Compton tail. M33 X-8 was well fit by this 
model with a Ax Vdof =1659.9/1534(1.08). Thus, these 
sources are still consistent with the black hole accretion 
model. Replacing the tbabs model used for the galac- 
tic column density with the tbvarabs, we measured the 
column density, oxygen abundance, and iron abundance. 
With this thermal disk and Comptonization model we 
obtained n^ = x 10^^ cm"^, 0/H = IM^oit 

and Fc/H = 2MtoVe, with a AxVdof = 1577.3/1534 
(1.03). Within the error bars, these results are consis- 
tent with those seen in Table El 

5. PROPERTIES OF THE ISM IN ULX HOST GALAXIES 

The major question to be examined in using ULXs as 
probes of the ISM is whether the hydrogen column den- 
sity and element abundances are primarily from the host 
galaxy or intrinsic to the local environment of the ULX. 
Before we can answer this question, it is important to 
understand the intrinsic spectrum of the source. In the 
previous section we discussed the nature of the soft com- 
ponent in light of the high signal-to-noise spectra of the 
14 ULX sources we examined. We found that if the spec- 
trum is due to thermal emission from a disk, modeling the 
spectrum with a variety of disk models (grad, diskbb, 
diskpn, bbody) does not significantly change the mea- 
sured oxygen abundance or hydrogen column density. 

Assuming the reliability of the hydrogen column den- 
sity and abundance measurements, based on their model 
independent values, we investigate the source of the ab- 
sorption in ULX spectra. In order to determine whether 
the model nn values suggest the necessity of extra local 
absorption, we compare the model values with column 
densities obtained from H I studies. We investigate this 
in section 5.1. The oxygen abundances (as an indica- 
tion of mctallicity), which we examine in section 5.2, can 
provide further clues of whether the absorption we see 
in the X-ray spectrum is intrinsic to the source. We also 
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Fig. 5. — Hydrogen column density from Eb—v (triangle) or H I 
studies (circle) vs. the hydrogen column density obtained from the 
tbvarabs model. H I column densities were obtained for Holmberg 
II, NGC 5204, and NGC 4559 through the WHISP survey (Swaters 
et al. 2002). Additional H I column densities were obtained for 
NGC 247, Holmberg IX, and M81 through VLA data (Braun 1995). 
The X-ray spectral fit columns are not biased towards significantly 
higher values than the alternate method column densities. This 
implies that most of the matter in the line of sight is from H I. 
This is not true, however, for the source M81 XMMl (represented 
by the two outlying points in the lower left portion of the graph) 
where the X-ray column density is much greater than those from 
the optical and radio. 



examine possible connections between the host galaxy's 
star formation rate and elemental abundances. 

5.1. Column Densities 

To determine whether the ULX X-ray hydrogen col- 
umn densities represent largely galactic column densi- 
ties or column densities local to the ULX, we com- 
pared the X-ray values to those obtained from opti- 
cal and radio studies. For a comparison to hydrogen 
column densities from optical studies, we used inter- 
stellar reddening values. In a study of dust scattering 
X-ray halos surrounding point sources and supernova 
remnants, Prcdehl & Schmitt (1995) derived a relation- 
ship between hydrogen column density and interstellar 
reddening, using X-ray data. They found that nn — 
5.3 X 10"^^ cm^^EB-v- They also found that these X-ray 
derived column densities are not affected by the intrinsic 
absorption of the X-ray source. Thus, the optical redden- 
ing, Eb-Vi becomes a useful tool in checking our own 
X-ray derived column densities. Through a literature 
searc h we found Eb_v values for th e sources M33 X-8 
(0.22:lLong. Charles, fc D^Ib^ l|20f)l ) and M81 XMMl 
(0.23; ;Kong_et al. (200(3)). The corresponding Eb-v- 
derived njj values are plotted in Figure [S] as triangles. 

For a comparison of X-ray derived hydrogen column 
densities with radio values, we obtained H I column den- 
sities for four objects (Holmberg II XMMl, NGC 4559 
X7, NGC 4559 XIO, and NGC 5204 XMMl) from the 
WHISP catalog l|Swaters et ai!.ll2002D . These are radio 
H I column densities within the host galaxy (galactic), 
not the Galactic/Milky Way columns. Exact values of 
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TABLE 6 

Column Densities from Optical 
AND H I Studies 



Source 




Method'' 


NGC 247 XMMl 


0.316 


H I 


M33 X-8 


0.170 


Eb-v 


Holmll XMMl 


0.157 


H I 


M81 XMMl 


0.122 


Eb-v 


M81 XMMl 


0.078 


H I 


HolmIX XMMl 


0.148 


H I 


NGC 4559 X7 


0.357 


H I 


NGC 4559 XIO 


0.300 


H I 


NGC 5204 XMMl 


0.182 


H I 



^Hydrogen column density (galactic, 
not Milky Way) determined by the 
method indicated in units of 10^^ cm~^ . 

Column densities computed from red- 
dening values (Eb—v) or from radio H 
I measurements, see section 5.1 for de- 
tails. 



the H I column densities were computed and given to us 
by Rob S waters. Additionally, we include H I column 
densities of NGC 247 XMMl, M81 XMMl, and Holm- 
berg IX XMMl from iBrauJ l|1995|) . We obtained the 
FITS files of H I column density maps from this paper 
(available on NED), where the pixel value corresponds to 
the galactic column density in units of lO^^cm^^. The 
H I n// values from both studies are plotted in Figure |5l 
as circles. The column densities, from the radio and red- 
dening studies, are listed in Table El 

We find that the host galaxy column densities from al- 
ternate methods (optical or radio studies) are not signifi- 
cantly different from the X-ray column densities. Partic- 
ularly, the X-ray values are not skewed towards substan- 
tially higher values than the optical/radio values. Thus, 
the X-ray columns are likely the galactic values without 
any additional local absorption. The exception, how- 
ever, is M81 XMMl (represented by 2 points in Figure[Sl 
one for each of the methods). The Eb-v value (1.22 x 
10^^ cm~^) and the H I value (0.78 x 10^^ cm~^) are sig- 
nificantly lower than the X-ray column density. This 
may indicate the presence of extra absorption around 
this source. 

The result that the X-ray hydrogen column densities 
are in good agreement with those from H I studies and in- 
terstellar extinction values is interesting considering that 
the X-ray measured column densities are along a direct 
line of sight to the ULXs while the H I measurements 
are an average over a larger beam area. The agreement 
between the two measurements implies that the ULX 
sources, with the exception of M81 XMMl, lie within 
roughly normal areas of their host galaxies (i.e. not 
in regions of higher column density such as a molecu- 
lar cloud). 

5.2. Elemental Abundances 

5.2.1. Test for Oxygen Ionization Level 

Before discussing implications of the determined oxy- 
gen and iron abundances, we relate a test performed to 
determine whether we could distinguish between differ- 
ent ionization levels of oxygen. To do this, we used 



the absorption edge model, edge, in XSPEC (using 
the full model: tbabs*tbvarabs*edge*edge* (grad + 
pow) ). We first checked to see that the abundance values 
obtained with the edge model matched the values from 
the tbvarabs model. We fixed the oxygen abundance in 
the tbvarabs model at zero and added the edge model, 
allowing the threshold energy and absorption depth (r) 
to float as free parameters (with an initial energy set to 
0.543 keV). We fit this model to the longest observations 
of Holmberg II XMMl and Holmberg IX XMMl in ad- 
dition to a source with a lower number of counts, NGC 
5408 XMMl (0112290601). 

To find the oxygen column density (no), we used the 
relationship that r = tr x no , where a is the cross section 
for photoabsorption. We used the cross section values for 
neutral oxygen published in Rcil man & Mans on (1979i) 
as an estimate. This choice is suppor ted by the results of 
lJuett. Schulz. fc Chakrabartvl l)2004f) . who measured the 
ratio of oxygen ionization states in the ISM as O II/O If» 
0.1. In Table the hydrogen column density, threshold 
energy, and optical dept h are listed for the sources fit 
with this model. From iReilman fc MansonI l)1979j) we 
used the cross section values of cr = 5.158 x 10~^^cm^ 
(for E = 0.540 keV) and a ^ 4.789 x lO^i^cm^ (for 
E — 0.570 keV). As seen in Table |3 the hydrogen column 
densities and the oxygen abundances obtained from this 
model are close to those from the tbvarabs model. The 
[0/H] values, where [0/H] = 12 + log(0/H) and O and 
H represent oxygen column density and hydrogen column 
density respectively, between the two models vary by less 
than 1%. 

To test whether the threshold energy from the edge 
model is affected by the ionization level of oxygen, we 
simulated spectra of an absorbed power law model with 
an oxygen edge. Tim Kallmann (P.C.) provided us with 
an oxygen ed ge model incorporating the cross sections of 
iGarcia et ali (j2005). The model allows for a variation of 
the ratio of O II/O I. Using the response and ancillary 
response matrices from the long Holmberg II ULX ob- 
servation, we simulated spectra with the XSPEC f akeit 
command for a F = 2.35 power law. We simulated spec- 
tra for O II/O I ratios of 0.0, 0.2, 0.4, 0.6, 0.8, and 
1.0. The simulated spectra were binned with a mini- 
mum of 20 counts/bin. The absorption edge component 
of the spectrum was then fit with the edge model. The 
fits to the threshold energies for the simulated spectra 
yielded values ranging from 0.53-0.59 keV, with no pref- 
erence of lower O II/O I ratios corresponding to lower 
threshold energies. Since the O I absorption edge occurs 
at E= 0.543 keV and the O II absorption edge occurs 
at E= 0.57 keV, we could not distinguish between these 
ionization states of oxygen. Our simulations show that 
the oxygen edge measurements will be sensitive to O I 
and O II but not to high ionization states (for instance 
O VIII which has an edge energy of 0.87keV). 

5.2.2. X-ray/optical [0/H] Comparison 

As noted above, we tested the oxygen abundances ob- 
tained with the absorption model tbvarabs against the 
abundances obtained from adding a photo-electric ab- 
sorption edge model, for three of the ULX sources. We 
found good agreement (< 1% difference in [O/H] values) 
between both models for the X-ray spectra. However, 
we found that it is not possible to distinguish between 
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Fig. 6. — Oxygen abundances from studies of H II regions vs. 
oxygen abundances from the tbvarabs model. The [O/H] ratios 
are given for Holmberg II XMMl, NGC 5408 XMMl, Holmberg IX 
XMMl, NGC 4559 XMMl, NGC 4559 XMM2, NGC 1313 XMM3, 
M33 X-8, NGC 300 XMMl, NGC 253 XMM2, and M81 XMMl 
using the P-calibration method (circle). Also, [O/H] ratios are 
included for M33 X-8, NGC 300 XMMl, M81 XMMl, NGC 253 
XMM2, and M83 XMMl using the R23 calibration method (tri- 
angle). Arrows indicate that the lower limit for the [O/H] X-ray 
parameter is below the plotted graph region. Our values are high 
compared to the P-method but in good agreement with the R23 
method. The sources with the largest difference between [O/H] 
values are those located in irregular galaxies. 



different low ionization states of oxygen using the edge 
model. 

We now discuss comparisons of our X-ray oxygen ab- 
sorption values with measurements in different wave- 
lengths, based on a hterature search for [O/H] ratios. 
In Figure El we compare our [O/H] ratios with those of a 
study conducted bv lPilvuein. Vzlchez. fc C ontinil (|2004D 
(circles). They provide a compilation of [O/H] ratios 
determined through spectrophotometric studies of H II 
regions. Their [O/H] values are based on the radial distri- 
bution of oxygen abundance using the P-method. They 
determined [O/H] values for spiral galaxies where pub- 
lished spectra were available for at least 4 H II regions. 
In addition, they reference [O/H] values for irregulars 
obtained through alternate methods. Our [O/H] val- 
ues are determined from the oxygen abundances listed 
in TableU Thus, [0/H]= 12 log(0 x 0.00049). O is 
the oxygen abundance obtained from the model, which is 
multiplied by the Wilms relative abundance of 4.9 x lO""* 
of oxygen to hydrogen. We were able to compare [O /H] 
values for the sources located in M33, NGC 253, NGC 
300, 1VI81, Holmberg II, NGC 4559, and NGC 5408. We 
include the [O/H] val ue computed for the Holmberg IX 
ULX bv llVrilleii ljl995|) of 8.12. This value was computed 
from an optical study of the surrounding H II region. 
Also, we add the [O/H] value of 8.4 for NGC 13 1 3, de- 
ter mined separ ately by both iCalzetti fc KinnevI l)1994f) 
and lWalsh fc"R ,ov (199'$. 

As shown in Figure El our [O/H] values are 
consistently high compared to those obtained from 
the iPilvugin. Vzlchez. fc Continil l)2004|) H II study. 



iPilvugin. Vzlchez. fc Continil l)2004|) include a discus- 
sion of how their values, obtained by the P-calibration 
me thod, are significantly lower than those obtained by 
jCarnett (2002) using the R23-calibration method. In Fig- 
ure we include [O/H] val ues for N GC 253, NGC 300, 
M33, 1VI81, and IVI83 from iGarnettI Ip0021 (triangles). 
Our oxygen abundances are in much better agreement 
with the values from this R23-calibration method. 

Our [O/H] values, which are consistent in the 
X-ray band between two separate absorption mod- 
els, (tbabs*tbvarabs*edge* (grad + pow) and an 
tbabs*edge*edge* (grad + pow) model), are consis- 
tent with the values of IGarnettI l|2002l) . We fur- 
ther wish to compare them to the metallicity pre- 
dicted from Sloan Digital Sky Sur vey (SPSS) r esults . 
As a result of a SDSS study, .Tremonti et o.l\ (|2004D 
found a luminosity-metallicity relation for their sam- 
ple of star- forming galaxies of: 12 -I- log{0/H) = 
-0.185(±0.001)Afs + 5.238(±0.018). In Figure (3 we 
compare our values with their results (represented by 
the line). We obtained absolute magnitudes (Mb) for 
the host galaxies using the total apparent corrected B- 
magnitude recorded in the HyperLeda galaxy catalogue 
(Paturcl et al. 198 ^) (parame ter btc) and the distances 
hsted in Table 1 of iPaoer ll For M33 and NGC 4559, 
which were not included in the previous study , we used 
the d istances of O.TlVIpc and 9.7 Mpc from iHo et al] 
l)1997|) . The majority of our sources are consistent with 
the SDSS results. However, the sources in irregular 
galaxies have metallicities much higher than predicted. 

5.2.3. Galaxy Properties 

More luminous galaxies are sometimes expected to 
have higher star formation rates, and thus higher metal- 
licity. However, we found no evidence of this. Investigat- 
ing further into the relationship between star formation 
rate (SFR) and metallicity, we chose to look at a galactic 
luminosity diagnostic that is less dependent on extinction 
from du st, the infrared galactic luminosity (hpin). In 
iPaper l'. we calculated L pj using data from the Infrared 
Astro nomical Satellite and the approach of Swartz et al\ 
1)20041]) . We quoted values of hpiR in fPaper ll Using the 
same method, with IRAS fluxes obtained from lHo et all 
|T99l . we find 1VI33 to have LpiR « 1.7 x lO^^ergs"! 
and NGC 4559 to have LpiR. « 7.4 x lO^^gj-gg-i^ 
used these L pi r values to compare the SFR to both the 
oxygen and iron abundances. 

In FigurejHl we plot oxygen abundance relative to SFR. 
All sources with the exception of Holmberg IX XMMl, 
which did not have available IRAS data, are plotted. We 
see, as was also illustrated in Figure that the lumi- 
nosity of the host galaxy does not determine the metal- 
licity. The more luminous galaxies do not have metal- 
licities higher than the less luminous galaxies. In fact, 
most sources have oxygen abundances that are roughly 
the Wilms solar abundances (indicated by the line). 

In Figure El we see that there is no relationship be- 
tween iron abundance (see values in Tables El and El for 
values) and SFR. The Wilms solar abundance for iron is 
only 2.69 x 10~^ relative the hydrogen abundance. The 
plots show that the metallicity relationship is very flat, 
all of the sources have roughly solar abundances. This 
is also seen in Figure ^1 This plot shows the [Fc/H] 
ratios versus the [O/H] values, both obtained through 
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TABLE 7 

Spectral Fits for Oxygen Edge with tbabs*tbvarabs*edge*edge*(grad + pow) model 



Source 




Eb 






[0/H]<= 


xVdof 


counts' 


Holmll XMMl 


l6+° oi 






7.2 


8.65/8.64 


2543/2017 


342874 


HolmIX XMMl 




0.543«:«1 


62+"" 


12.0 


8.81/8.82 


2887/2751 


148061 


NGC 5408 XMMl 


-0.03 


0.538+0Q2 


Q 07 + 0.29 

'^■•^'-0.21 


7.1 


8.61/8.63 


298/334 


10045 



^Hydrogen column d ensity determined fr o m tbv arabs in units of lO'^^ cm ^. The Galactic value 
of n^f was fixed to the lDickev fc: LockmanI 1119901) value with the tbabs model. 

''Threshold Energy obtained from the edge model in keV. Note that one edge model was used to 
correct for the difference in the oxygen edge between the pn, MOSl, and MOS2 CCDs. The other 
edge model was used to measure the oxygen abundance from the 542 eV K-shell edge. 
Absorption depth obtained from the edge model 

■^Column density of oxygen estimated from the edge model in units of 10^^ cm~^ 

Abundance of oxygen relative to hydrogen from the edge model versus the value quoted in Tablel2l 
[O/H] = 12 + log(0/H) where O is oxygen abundance and H is hydrogen abundance 

'Total number of photon counts from the pn and MOS detectors 




Fig. 7.— Host galaxy Mg vs. ULX X-ray [O/H] ratio. The 
solid line represents the SDSS results from a study of star-forming 
galaxies (.Tremonti et al. 2 003) . Arrows indicate that the lower 
limit for the [O/H] X-ray parameter is below the plotted graph 
region. Our values are largely consistent with the Sloan results, 
with the exception of Holmberg IX and Holmberg H. Both of these 
galaxies are irregulars. The X-ray [O/H] ratio derived from the 
ULX spectra is higher than predicted for these objects. 



Fig. 8. — Oxygen abundance ([O/H]) as a function of the host 
galaxy's FIR luminosity, obtained from IRAS. The line represents 
the Wilms solar value of [O/H] . Arrows indicate that the lower 
limit for the [O/H] X-ray parameter is below the plotted graph 
region. There appears to be no correlation between hpfn and 
oxygen abundance. In fact, the sources show roughly solar [O/H] 
abundances regardless of host galaxy. 



the tbvarabs modeL The solar Wilms values are repre- 
sented on the plot by the open circle symbol. It appears 
that the sources are slightly more abundant in iron than 
the solar value, however, the error bars are quite large. 
The oxygen abundances, as before stated, are roughly 
the solar Wilms value. Roughly, the abundances appear 
solar. 

With such a flat relationship between abundance and 
luminosity, we tested to see if this result carried through 
in a comparison of abundance versus radial distance 
within the host galaxy. This follows upon an interesting 
property observed in spiral galaxies. Namely, that abun- 
dances are typically higher in the ce nter of the g alaxy 
and decrease with increasing radius llSearlel 119711) . We 



tested how our results compare to this result in Figure 
111! For the sources located within spiral galaxies, we 
plotted the [O /H] ratio as a function of distance R from 
the dynamical center (as repo rted by N ED). We used the 
distances quoted in Table 1 of lPaperTI to translate angu- 
lar distance on the sky into kpc from the galactic center. 
As evidenced in the plot, we do not see much variation in 
the [O/H] values. Clearly, the expected scaling of higher 
oxygen abundance towards the center is not seen. Since 
the sources lie in different host galaxies with different 
relative abundances, it is possible that this trend may 
be detectable with a larger sample of galaxies. However, 
the implication from Figures □ through ini is that the 
environment of the ULX sources is relatively uniform in 
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Fig. 9. — Iron abundance ([Fc/H]) as a function of the host 
galaxy's FIR luminosity obtained from IRAS. Holm IX did not 
have corresponding IR data, however, being an irregular galaxy, its 
hpfn would probably be comparable to Holm II. The iron abun- 
dances of both of the irregulars are well above those of the ULXs 
in spiral galaxies. The sources in galaxies with larger L^/jj have 
roughly equal iron abundances, within the error bars. The iron 
abundances in all cases are well above the Wilms solar abundance. 
Arrows indicate that the lower limit for the [Fe/H] X-ray parameter 
is below the plotted graph region. 
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Fic;. 11. — Oxygen abundances from the tbvarabs model vs. 
distance of the ULX from the host galaxy's dynamical center. 
The [O/H] ratios are given for spiral galaxies only. There is no 
correlation between radius and oxygen abundance in our sample. 
However, since the sources originate in different host galaxies with 
different abundances, there is no conclusion that can be drawn 
from this data. The host galaxy NGC 4559, containing two ULX 
sources studied, does appear to show a possible correlation (the 
source nearer the center has a higher oxygen abundance) , however, 
the error bars on the oxygen abundance are large for both sources. 
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Fig. 10. — [Fe/H] as a function of [O/H] from the tbvarabs 
model, using the grad and power law fit. The Wilms solar abun- 
dance is indicated with an open circle. The ratio of Fe/O abun- 
dances obtained through the X-ray spectral fits are approximately 
the Wilms solar values. Arrows indicate that the lower limit for 
the [Fe/H] X-ray parameter is below the plotted graph region. 



terms of metallicity. The ULX sources appear to live in 
similar environments, with metallicities roughly solar. 

6. SUMMARY 

Through our work, we conclude that X-ray spectral 

fits to ULX sources do provide a viable method of find- 
ing abundances in other galaxies. We have determined 
hydrogen column densities and oxygen abundances along 



the line of sight to 14 ULX sources. To measure these val- 
ues, we assumed a connection between ULXs and Galac- 
tic black hole systems such that the ULX spectra used 
in this study correspond to a high/soft state. Therefore, 
we modeled the sources with an absorbed disk model and 
power law (in XSPEC tbabs *thva.Ta.hs* (accretion disk 
model + pow)). We tested the effect different accretion 
disk models have on the measurement of the host galaxy's 
hydrogen column density and elemental abundance with 
the disk models: grad, diskpn, diskbb, and bbody. We 
found that the measured hydrogen column density and 
abundances are model independent. 

We also tested the physical plausibility of this model by 
comparing the mass and mass accretion rates obtained 
from the grad model with expected results based on 
Galactic black hole systems. We found that the ULX 
spectra were consistent with the high/soft state, with 
L/LeiIcL values < 1.0 for sources with a standard ther- 
mal model at low energies and power law dominated 
higher energy spectrum. For four sources, the spectra 
were consistent with a heavily Comptonized spectrum. 
These sources are more likely stellar mass black hole sys- 
tems in a very high state of accretion. We modeled their 
spectra with a power law at low energy and disk model 
around 1 keV. This model provided similar column den- 
sity and abundance values to a more physical absorbed 
disk and Comptonization model. 

Comparing our X-ray measured column densities with 
those from optical and H I studies for 8 of our sources, 
we find that 7 of the sources have X-ray column densities 
approximately equal to those of the alternate methods. 
This implies that the hydrogen columns towards most 
of our ULX sources represent that of their host galaxy. 
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Since the H I studies represent averages over a large beam 
area where the X-ray column densities are directly along 
the line of sight to the ULX source, this implies that the 
ULX sources lie within roughly normal areas of their host 
galaxies. The exception in this study was M81 XMMl, 
whose X-ray hydrogen column density was large relative 
to the H I study. This suggests that there is extra ab- 
sorption intrinsic to this source. 

The oxygen abundances appear to be roughly the 
Wilms solar values. For five sources, the count rates were 
sufficient to determine iron abundances without large er- 
ror bars (see Figures El and EU- We found that iron 
abundances for these sources were slightly overabundant 
relative to the solar Wilms value. However, within error- 
bars, the abundances appear solar. X-ray derived [0/H] 
values are com parable to those from an optical study by 
l)Garnettll2002() . indicating that the X-ray derived values 
are the same as the [O /H] values of H II regions within 
the host galaxy. Luminosity-metallicity relationships for 
the ULX host galaxies show a flat distribution, as does 



a radius-metallicity plot. Therefore, it appears that the 
ULX sources exist in similar environments within their 
host galaxy, despite the wide range of host galaxy prop- 
erties. 
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Fig. A12. — The number of counts for simulated spectra versus the [O/H] value (left) and [Fo/H] value (right) from the 
tbabs*tbvarabs* (grad + pow) , where errors represent the 90% confidence rate. The horizontal lines represent the solar Wilms values 
([0/H]f5i 8.69 and [Fe/H]Ri 7.43). Arrows are used to represent errors that ex tend below the range of the plot. The tip of the arrow point 
represents the [O/H] or [Fe/H] parameter from the XSPEC model (see Table IMl. 



APPENDIX 
SPECTRAL SIMULATIONS 

We conducted spectral simulations in order to: (1) determine the number of counts needed to measure the oxygen 
and iron abundances and (2) to verify the model independence of the galactic column density and abundances with 
respect to the grad, diskbb, and bbody models (seen in a comparison of Tables |21 and 2J). Towards this end, we created 
simulated pn spectra based on the long Holmberg IX XMMl observation's (0200980101) unbinned, pn spectrum. We 
used the base (grad + pow ) model parameters as indicated in Tabled We modeled both the Galactic column density 
ijDickev fc Lockman (1990) value: hh — 4.0 x lG^°cm~^) and host galaxy column density [hh = 1-9 x 10^^ cm^^) 
with individual tbabs models. Thus, all of the abundances were set to the solar Wilms values. We used the XSPEC 
command f akeit to create simulated spectra with 200000, 40000, 10000, 5000, and 2000 counts. The simulated spectra 
were binned with 20 cts/bin using grppha. We fit the binned simulated spectra with the models tbabs*tbvarabs* (grad 
+ pow), tbabs*tbvarabs* (diskbb + pow), and tbabs*tbvarabs* (bbody + pow). This allowed us to see the effects 
the different models have on the measured galactic hydrogen column density and abundances. The results for these 
fits are seen in Table IXSl The range that the oxygen and iron abundances were allowed to vary within was 0.0 (lower 
limit) to 5. (up per limit) with respect to the solar values. 

In Figurc [Al2l the number of simulated counts versus the errors on the oxygen and iron abundances are plotted for the 
tbabs*tbvarabs*(grad + pow) model. Here, [0/H]= 12 + log(0 x 0.00049) and [Fe/H]= 12 + log(Fe x 0.0000269), 
using the solar Wilms values for O/H and Fe/H. As seen in the plots, the errors in oxygen abundance are much 
smaller for a given number of counts compared to the errors in iron abundance. Further, the upper limits on the 
oxygen abundance continue to be meaningful through 2000 counts. This is not true for the iron abundances, where 
the error bars extend through the entire range of allowed values (from Fe/H = 0.0 - 5.0, or [Fe/H] up to 8.13). 
Thus, our simulations show us that the iron abundance (from measurements of the Fe L-shell edge at 851 eV with the 
tbvarabs model) requires at least 5000 counts to be detected. At 5000 counts the model derived value (Fe/H — 1.36) is 
meaningful, however, the errors extend throughout the entire allowable range (Fc/H ~ 0.0-5.0). The oxygen abundance 
(from measurements of the O K-shell edge at 542 eV with the tbvarabs model) is detected down to 2000 counts, but 
with large err ors b elow 10000 counts. 

From Table 1X51 we find that the same trends described in Section 4.1 are present in our simulations. Namely, there 
is little variation between the model derived abundances and column densities. A comparison of the mean and 
oxygen abundance values shows « 2% difference between the grad and diskbb model. The bbody nn values are 
roughly 26% higher while the oxygen abundances are w 11% higher. Comparing the 40000 and 200000 count spectra 
for the iron abundance, we find that the grad and diskbb model values differ by « 12% while the bbody model results 
are larger by a factor of 50%. While the bbody model yields lower column densities and higher abundances, the diskbb 
and grad models are in agreement. The differences in the bbody results are low for the column density and oxygen 
abundance, but appeared significant for the iron abundance. 
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TABLE A8 

Spectral Fits to Simulated HolmIX XMMl tbabs*tbabs* (grad + pow) 

Spectrum 



Counts'' II jf'' Oxygen abundance'^ Iron abundance'^ x^/<iof 



tbabs*tbvarabs* (grad + pow) 



200000 


0, 


91 +0.02 
■^^-0.02 


i.08ig:« 


07+0.68 
"■*''-0.76 


1687/1727 


40000 


0, 


-IQ+0.04 
■-'■'-0.03 


1 ic; + 0.26 
'■'3-0.33 


"■'"-0.97 


933/914 


10000 


0, 


9n+o.07 

■^^-0.06 




n Qn+2-64 
"■y"-o.9o 


340.5/401 


5000 


0, 


16+°" 
■-'"-0.07 


1 07+0.75 
'■•''-1. 34 


'■•J"-!. 36 


187.2/201 


2000 


0, 


nq+0.14 

■•'•'-0.19 


1 20+°-" 
'■^"-1.20 


"■^-5.0 


71.9/89 


tbabs*tbvarabs* (diskbb + pow) 


200000 


0, 


91 +0.02 
■^^^-0.02 


'■Ll»-0.19 


og+0.94 
"■'"'-0.76 


1691/1727 


40000 


0, 


1Q+0.03 
■'^-0.03 


'■"'-0.35 


1 oq+1-4^ 

'■""-1.09 


934/914 


10000 


0, 


1Q+0.03 
■'''-0.05 


1 cro + 0.43 
l-^°-0.50 


1 01+2-*° 
'■"'-1.01 


341/401 


5000 


0, 


■'"-0.08 


1 41+0.76 
'■^'-1.16 


1 c:7+3.43 
'■•"-1. 57 


187.5/201 


2000 


0, 


09+0. 18 

■•J^-o.io 


1 90+0.59 
'■^•'-1.23 


"■"-5.0 


71.7/89 








tbabs*tbvarabs* (bbody 


- + pow) 




200000 


0, 


,-+0.02 

■" -0.02 


'■^^-0.20 


'■""-1.01 


1702/1727 


40000 


0, 


-|c:+0.04 
■'•5-0.02 


1 90+0.36 


1 70+1.97 

'■ '•'-1.72 


937/914 


10000 


0, 


-|c:+0.07 
■'•'-0.04 


1 77+0.61 
'■ ' ' -0.72 


-1 99+3.78 
'■^^-1.22 


343.7/401 


5000 


0, 


1 0+0.10 
■'•'-0.07 


'■"'-i.ei 


9 1 Q+2.81 
z.l»_2 jg 


187.5/201 


2000 


0, 


91 +0.19 
■^'-0.12 


'■'"-1.16 


"■"-5.0 


71.8/89 



^Total number of photon counts for simulated pn spectrum. Simulated spectra 
created with the XSPEC f cikeit command, using the model parameters and response 
files from the long HolmIX XMMl observation. 

^Hydrogen column density determine d from tbvarabs in u n its of 10^^ cm~^. The 
Galactic value of nn was fixed to the iDickev 8z Lockma value with the 

tbabs model. 

^Element abundance relative to the Wilms solar abundance from the tbvarabs 
model 
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TABLE B9 

Spectral Fits for EPIC spectra with tbabs*tbvarabs*(grad + pow) model 



Source 




O abund.'' 


Fe abund.'' 


Mass (M©) 


M/MEdd" 


F 


xVdof 


counts'^ 


HoImIX XMMl 


U-lO-0.02 


-■-•0' -0.14 


1-D1_0.88 


442tJiO 




1 47+0.02 

l-*'-0.02 


1637.1/1532 


171210 




U--'-»_0.04 




1 S4+2'23 
-^■^^-1.84 






1 fiS+0-08 
1-O»_0.09 


477.8/486 


14996 


Holmll XMMl 


12+'''0i 


-, 10+0.14 
-^•-^•^-0.14 


0+0-24 




11+0-03 
'-'•^^-0.03 


9 45+0.06 
^•^^-0.05 


939.3/932 


150290 






1 iq+0.51 


3-08^3^08 


^'^-163 


'-'•'^'-0.03 


9 01+0.12 
^••^1-0.13 


362.4/368 


13983 


M33 X-8 


o.oslH! 


1 91+0.37 

-'-•^-'--0.48 


'-'■^-0.0 


7 90+0.68 
'•^'-'-0.63 


'-'•^"-0.06 


9 19+0.25 
^•1^-0.26 


714.9/749 


46579 


M81 XMMl 


0.48l»;« 


1 99+0.11 
-■-•^^-0.15 


1 61+°'5* 
1-O1-1.09 




90+0.18 
■^•^^-0.31 


01 +0.68 
■^•■^1-0.95 


865/811 


42188 


N4559 X-7 


16+'^-°'' 
'^■^°-0.03 


'J-0^-0.52 


OQ+1.08 


lOOOt^^f 


'-'•'^^-0.03 


^•l"-0.08 


407.1/359 


12492 


NGC5204 XMMl 


'^■"^-0.01 


l-28ti'.28 


"■"-0.0 


,90+146 

l^U_gQ 


04+0-05 
'-'•'^^-0.02 


1 04+0-10 
l-y^-0.13 


305.7/323 


9991 




09+° °3 
U-U»_o.o3 


^•"'^-0.99 


0+2-*^ 
'J-'^-o.o 






l.MU_Q 25 


244.8/278 


7503 


NGC4559 X-10 


n n+"-03 

'J-ll-0.03 


1 47+0.63 
-"-•^'-0.53 


"■"-0.0 




4.76t^:i 


9 99+0.58 
^•^^-0.29 


266.4/328 


8834 


NGC300 XMMl 


09+'^ °=* 
'^■'-'^-0.03 


9 co+0.81 


OQ+1.67 
"■"-0.0 


leotgf 


^•'-'^-0.01 


2 46+°-^5 
^•^'-'-o.ie 


278.9/237 


6771 


NGC1313 XMM3 


40+"-^'^ 




o r, + 1.46 

■'■'J-1.77 


907+1477 
^'^'-695 


^•'^"-o.oi 


9 9c:+0.09 
^•^^-0.17 


246.2/232 


6237 


NGC5408 XMMl 


0.05tH^3 


2.43ti:g 


^■'J-s.o 




11+0.04 
'-'•11-0.03 


9 91+0.31 
^•^^-0.30 


174/187 


5928 


NGC4631 XMMl 


0.28l»:- 


0.41«;- 


'J-^-o.o 




72+1.15 
'J- '^-0.72 


1 05+0-24 
l-»J„0.40 


215.4/189 


4830 



^Hydrogen colum n density determined from tbvarabs in units of 10^^ cm ^. The Galactic value of n/j was fixed to the 
iDickev & LockmanI il990f) value with the tbabs model. 
^Element abundance relative to the Wilms solar abundance from the tbvarabs model 
^Ratio of mass accretion rate from the grad model to Eddington accretion rate (see Section 4) 
*^Total number of photon counts from the pn detector 



COMPARISON WITH SAS 6.5 DATA PROCESSING 

We reprocessed the pn data in this study with the new version of SAS (6.5), in order to observe the effects different 
SAS versions would have on our results. After processing the data (with epchain), we extracted spectra as described 
in Section 2. The spectra were then fit with the tbabs*tbvarabs* (grad + pow) model. Note that the edge model 
was not needed to correct for differences in calibration between the MOS and pn since we only fit pn spectra. Since 
the MOS data was not added, the counts for the pn only spectra are lower than the MOS + pn sp ectra . Thus, we 
excluded spectra for sources with less than 5000 counts in the pn. The results are shown in Table IB9I ordered by 
number of counts. 

In Figure IIjl3l we plot a comparison of the X-ray hydrogen column densities and oxygen abundances (both from 
the tbvarabs model) from SAS version 6.0 and 6.5. From these plots, we find that the SAS 6.0 n^f values are slightly 
lower but consistent with the SAS 6.5 observations. The [O/H] values are slightly higher for SAS 6.0, but are also 
consistent. Since the hydrogen column density and oxygen abundance values are consistent between both versions, 
we did not find it necessary to recalculate EPIC spectral parameters for the sources. Thus, the values quoted in the 
paper are from the data processed with SAS version 6.0. While the error bars for the SAS 6.5 Fe abundances are 
larger, due to the lower number of counts in the pn alone, within the errors the values are consistent with the SAS 
6.0 values. For the highest number of count sources (Holmbcrg IX XMMl and Holmberg II XMMl), the fits for the 
longest observations yielded the same Fe abundance for Holmberg II XMMl (Fe/H = 0.0 ranging to « 0.20) and 
consistent values for Holmberg IX XMMl (the SAS 6.0 values (Fe/H from 2.19 to 1.94) were well within the SAS 6.5 
values (Fe/H from 2.37 to 0.75)). 

We also note that the release notes for SAS version 6.0 (Carlos Gabriel & Eduardo Ojero, available on- 
lin e) indicate no substantial d ifferences in the EPIC responses. Therefore, the edge model correction quoted 
in iBaumgartner fc Mushotzky| ( ^2006') is adequate and did not need to be recalculated from SAS 5.4 (as in 
IBaumgartner fc Mushotzkvl l|2006h ) to SAS 6.0. 
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Fig. B13. — Comparison of the galactic/ULX hydrogen colu mn density obtained from the tbvarabs model for SAS 6.0 processed data 
(from Table|3J versus SAS 6.5 processed data (from Table lB9l (left) and comparison of the similarly obtained [O/H] values (right). From 
these plots, we find that the SAS 6.0 nn and [O/H] values are consistent with the SAS 6.5 values. The version 6.0 column densities are 
slightly lower than the 6.5 values and the version 6.0 oxygen abundances are slightly higher, on average. The lines in the plot represent 
agreement between the SAS 6.0 and 6.5 values. Arrows are used to indicate fits where the lower limit for the oxygen abundance extends 
below the plotted range. 



